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ABSTRACT 

We present the kinematic properties of a tidally disrupted dwarf galaxy in 
the Milky Way, based on the hypothesis that its central part once contained the 
most massive Galactic globular cluster, uj Cen. Dynamical evolution of a self- 
gravitating progenitor galaxy that follows the present-day and likely past orbits 
of uo Cen is calculated numerically and the kinematic nature of their tidal debris 
is analyzed, combined with randomly generated stars comprising spheroidal halo 
and flat disk components. We show that the retrograde rotation of the debris 
stars at ~ —100 km s~^ accords with a recently discovered, large radial velocity 
stream at ~ 300 km s~^ towards the Galactic longitude of ~ 270°. These stars 
also contribute, only in part, to a reported retrograde motion of the outer halo 
at the North Galactic Pole. The prospects for future debris searches and the 
implications for the early evolution of the Galaxy are briefly presented. 

Subject headings: Galaxy: formation — globular clusters: individual {uo Cen- 
tauri) — stars: kinematics 



1. INTRODUCTION 

Omega Centauri, the most massive globular cluster in the Milky Way, is unique in terms 
of its metallicity content, internal kinematics, and structure. Unlike other Galactic globular 
clusters, uj Cen shows a wide spread in metallicity (e.g. Norris, Freeman, & Mighell 1996), 
with a main metal-poor component at [Fe/H]~ —1.6, a second smaller peak at [Fe/H]~ —1.2, 
and a long tail extending up to [Fe/H]~ —0.5. The metal-rich population holds a low 
velocity dispersion and no sign of rotation, in contrast to the rotating metal-poor population. 
Furthermore, the metal-rich stars in uj Cen are largely enhanced in s-process elements relative 
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to those in globular clusters and field stars with similar metaUicities (e.g. Norris & Da Costa 
1995), thereby suggesting that the ejecta from low- mass, asymptotic giant branch (AGB) 
stars had to be retained and incorporated into the next-generation stars. 

In spite of its large mass (5 x 10^ M©), it has been demonstrated by Gnedin et al. (2002) 
that Lo Cen is not unique in its abihty to retain the AGB ejecta as found for other clusters. 
An isolated formation of u; Cen is thus unlikely, because the enriched gas would easily be 
lost by encountering the Galactic disk. The most viable explanation for the uniqueness of lo 
Cen is that it was once the dense nucleus of a dwarf galaxy (Freeman 1993). A gravitational 
potential provided by progenitor's stellar system and dark matter (as suggested from dwarf 
spheroidal galaxies in the Local Group, Mateo 1998) would help retaining the enriched gas 
and let the cluster being self-enriched at least over a few Gigayears. 

If this hypothesis is the case for the origin of uj Cen, the question arises: where and 
in what form does the stellar system of its progenitor galaxy remain? Dinescu (2002) first 
investigated this issue, by examining the possible signature of the progenitor's tidal debris 
among nearby metal-poor stars in the catalog of Beers et al. (2000, BOO). She identified a 
group of stars with —2.0 <[Fe/H]< —1.5, which departs from the characteristics of the inner 
Galactic halo but has retrograde orbits similar to cu Cen. Her simplified disruption model of 
the progenitor galaxy demonstrated that trailing tidal debris, having orbital characteristics 
similar to the chistcr, can be found in the solar neighborhood, although the concrete spatial 
distribution and kinematics of the debris stars remain yet unclear. 

This work motivates us to conduct an N-body simulation for the tidal disruption of cu 
Gen's progenitor galaxy, to obtain the characteristic structure and kinematics of its debris 

stars and compare with various observations showing signatures of recent merging events in 
the Milky Way (Gilmore, Wyse, & Norris 2002, GWN; Kinman et al. 2002, K02; Chiba 
& Beers 2000, CB). In particular, we show that a recently identified stream of stars at 
heliocentric radial velocity of ~ 300 km s^^ (GWN) is a natural outcome of the current 
disruption model, without affecting local halo kinematics near the Sun and microlensing 
optical depth towards the Large Magellanic Cloud (LMC). 

2. SIMPLE DYNAMICAL MODEL OF A PROGENITOR GALAXY 

We calculate the dynamical evolution of an orbiting dwarf galaxy in a fixed external 
gravitational potential representing the Milky Way. The potential consists of three parts: a 
spherical Hernquist bulge $6(r), a Miyamoto-Nagai disk ^d{R,z), and a logarithmic dark 
halo $/i(r), where r is the Galactocentric distance and (i?, z) denote cyhndrical coordinates. 
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Each is given as, ^^(r) = -GMb/{r + a), = -GMa/ \J + {h + ^/WT^f, and 

^f^{r) = vl/2\n{r'^ + (f), where = 3.4 x 10^° Mq, a = 0.7 kpc, = lO^^ M©, 6 = 6.5 
kpc, c = 0.26 kpc, t;^ = 186 km s~^, and d — 12 kpc. This choice yields a circular velocity 
of 228 km s~^ at the solar circle of it!© = 8 kpc and a flat rotation curve outside it!©. 

We set self-gravitating particles in the dwarf galaxy following a King model, where the 
central density, central velocity dispersion, and core radius are given as 0.3 M© pc~^, 18.1 
km s~^, and 0.56 kpc, respectively. In addition, a particle with the mass of 5 x 10^ M© 
representing uj Cen is placed at the center of the galaxy to trace its orbit. This setting yields 
the total mass of the galaxy as Mtot = 5.79 x 10^ Mq. A part of the mass is provided by stars, 
which is roughly estimated from the mean metallicity of stars in uj Cen (([Fe/H]) ~ —1-6), 
combined with the metallicity-luminosity relation for the Local Group dwarfs (Cote et al. 
2000) and the mass-to-hght ratio (assuming M/L ~ 4 obtained for u> Cen, Meylan et al. 
1995), yielding M^tars ~ 10^ M©. Thus, our model galaxy is largely dominated by a dark 
component, in agreement with the observed large M/L in dwarfs (Mateo 1998). The galaxy 
is represented by a collection of 10^ particles and the self-gravity is calculated in terms of a 
multiple expansion of the internal potential to fourth order (Zaritsky & White 1988). 

In the course of its orbital motion, a dwarf galaxy is disrupted by Galactic tides, whereas 
its dense core is expected to survive and follow cu Gen's orbit. While our calculation in a 
fixed Galactic potential neglects dynamical friction against progenitor's orbit, the effect 
is only modest for the system of < 10^ M© (Zhao 2002), especially dming a few orbital 
periods required for tidal disruption. We thus examine two representative orbits for the 
progenitor, model 1 and 2: model 1 follows the current orbit of uj Cen, whereas for model 2, 
we calculate an orbit back to the past over ~ 2 Gyr from its current position and velocity 
by fully taking into account dynamical friction and then set a progenitor galaxy on its non- 
decaying orbit. These two models provide us with satisfactory information on the generic 
properties of a tidally disrupted progenitor and we postulate that the realistic nature of 
their debris is midway between these model predictions. We calculate cu Gen's orbit, based 
on the distance^ from the Sun D = 5.3 ± 0.5 kpc, position (l,b) = (309°, 15°), proper 
motion (yU^ cos 5, /x^) = (— 5.08±0.35, — 3.57±0.34) masyr"^, and heliocentric radial velocity 
Vios = 232.5 ±0.7 km s~^ (Dinescu, Girard, & van Altena 1999). This orbit for model 1 
is characterized by frequent disk crossings with a period of Torb = 0.8 x 10^ yr, retrograde 
motion, and apo and pericentric distances {rapo,rperi) = (6.4,1.1) kpc. For model 2, we 
obtain Tgrb — 1-5 x 10^ yr and {rapo,rperi) = (11.3,3.0) kpc. In both experiments, we adopt 
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the same progenitor mass Mtot and place it at apocenter to maximize its survival chances . 

Figure 1 shows the spatial distribution of the tidally disrupted debris. Upper (middle) 
panel shows model 1 (model 2) after the 1.37 (1.86) Gyr orbital excursion of the progenitor 
galaxy. Lower panel shows the orbit of the galaxy center. In the course of the orbital motion 
of the galaxy, its structure is made elongated along the orbit induced by Galactic tides, in 
particular at its periccntcr passages, and then the particles are spread out to form the tidal 
streams along the orbit. A rosette-like feature of the debris becomes steady after about eight 
orbital periods. Model 1 results in more compact distribution than model 2, which reflects 
the difference in orbital radii. 

Figure 2 shows the velocity distributions of the debris particles in cylindrical coordinates 
(vr, v^, Vz). As is evident, model 1 and 2 provide essentially the same debris kinematics: most 
remarkable is a sharply peaked distribution at ~ —100 km s~^, arising from a retrograde 
orbit of a progenitor. These kinematics suggest that the difference in model 1 and 2 resides 
only in the spatial extent of the debris. 

3. EFFECTS OF DEBRIS ON GALACTIC KINEMATICS 

3.1. Method for kinematic analysis 

In order to assess the reality of the debris stars in light of observed stellar kinematics in 
the Milky Way, we analyze the kinematics of both the simulated debris and other Galactic 
stars generated randomly by a Monte Carlo method. The metal-poor halo is modeled as a 
flattened spheroid p oc (i?^-|-z^/g^)~^-^/^, where q is an axis ratio ranging 0.55-0.7, anisotropic 
velocity ellipsoid (ctr, cr,^, cr^) = (154,121,96) km s~^, and small mean rotation {v^) = 24 
km s~^, as found for halo stars with [Fe/H]< —2 near the Sun (CB). Thin and thick disks 
are modeled as p oc exp{—R/ R^i) sec'^ {z / z^) , where Rii = 3.5 kpc and z^ = 0.3 (1) kpc for thin 
(thick) disk. Both disks rotate at 200 km s~^, having velocity eUipsoids of (34, 25, 20) km s~^ 
and (46,50,35) km s~^ for thin and thick disks, respectively (CB). The relative fraction of 
each component is flxed using observed local densities near the Sun, in such a manner that 
the halo and thick-disk densities a,t D < 1 kpc are 0.2 % and 2 % of the thin-disk density, 
respectively (e.g. Yamagata & Yoshii 1992). 

In our model of cu Cen's progenitor galaxy, the self-gravitating particles represent both 
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stars and dark matter. We note that a correct estimate for the fraction of stars is uncertain, 
because their M/ L ratio as well as the amount of dark matter in the progenitor is unavailable. 
As a useful method to incorporate this ambiguity for the current kinematic analysis, we set 
a parameter / as the fraction of the debris particles relative to halo stars near the Sun, when 
all of the particles are regarded as stellar ones. By this, the normalization of the halo density 
is obtained for the given debris particles. 

A typical value of / for the conversion of the simulated particles to the stars is estimated 
in the following manner. Model 1 (model 2) yields 21 (74) particles aX D < 2 kpc, giving 
the mass density of pg = 0.4(1.3) x 10~^ M0 pc~^ near the Sun, whereas the total mass 
density and metal-poor halo density have been derived as 8 x 10~^ Mq pc^^ (Gates, Gyuk, 
& Turner 1995) and 6.4 x 10"^ Mq pc'^ (Gould, Flynn, & Bahcall 1998), respectively. 
Then, if the debris stars (with Mgtars ~ 10^ Mq) are distributed in the same manner as the 
simulated particles (with Mtot — 5.74 x 10^ Mq), which would be a reasonable approximation 
in view of the dissipationless nature of stars, the mass density of the debris stars in the solar 
neighborhood can be estimated as {M stars/ Mtot) Pg = 0(10^^) Mq pc^^, which is about 1 % 
of the halo density. Thus, /, defined here as debris fraction at D < 2 kpc, is expected to be 
of order of a few per cents. 



3.2. GWN's radial velocity survey 

Recently, GWN reported a spectroscopic survey of ~ 2000 F/G stars down to V = 19.5 
mag, in the direction against Galactic rotation (/, b) = (270°, —45°) and (270°, +33°), where 
radial velocities, vios, in combination with distances largely reflect orbital angular momen- 
tum. The Vios distribution of the stars a few kpc from the Sun shows two stellar streams at 
vios ~ 100 km s~^ and ~ 300 km s~^, which are not explained by known Galactic compo- 
nents. While the stream at vios ~ 100 km s~^ was reproduced by their model of a merging 
satellite in prograde rotation, the stream at vios ~ 300 km s""^ remains yet unexplained. 

Figure 3a shows the Vios distribution for the debris stars of model 2 and halo stars of 
q = 0.7 (i.e. without disks) at 1 < D < 5 kpc, 260° < / < 280°, and two fields for b. 
Figure 3b shows when disk stars are incorporated. As is evident, the debris stars from cu 
Cen's progenitor form a local peak at vios ~ 300 km s~^, which is provided by many stars 
having ~ —100 km s~^. This is in good agreement with the vios ~ 300 km s~^ stream 
discovered by GWN. A more flattened halo than the case g = 0.7 yields a higher peak, since 
the density contrast of the debris relative to halo is made higher in this survey region. It is 
worth noting that model 1 yields essentially the same w/os distribution as model 2, reflecting 
the same velocity distribution, although to attain the same peak height at ~ 300 km s~^, / 
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be a few factor larger and the selected range of Z be a few degree higher than the respective 
values in model 2, because of less number of debris stars near the Sun. This rule applies to 
other considerations below as well. 



3.3. Kinematics at the North Galactic Pole 

Majewski (1992) suggested that the outer halo at the North Galactic Pole (NGP) shows 
a retrograde rotation (i><^) ~ —55 km s""^ at 2; > 4 kpc. Also, K02 reported that their sample 
of horizontal branch stars at 2 < ^ < 12 kpc shows a retrograde rotation at {v^) ~ —65 
km s~^. On the other hand, halo stars near the Sun show no retrograde rotation (CB). 

To investigate the role of the debris stars in this issue, we select those of model 2 and 
randomly generated stars at 6 > 70° and 2 < D < 5 kpc (resembhng K02's selection). Since 
the observational determination of full space velocities involves rather inaccurate information 
of proper motions compared to radial velocities, we convolve the velocity distribution of stars 
with a Gaussian distribution for velocity errors, having 1 o" of a typical 30 km s^^ error. The 
resulting velocity distribution shows a non-Gaussian feature owing to the presence of the 
debris stars: the distribution holds an extra peak at ~ —100 km s^^ in addition to the 
~ 20 km s^^ peak, where the former amplitude becomes comparable to the latter one 
at / of a few percents, whereas for vr and Vz, the velocity distributions are made slightly 
asymmetric. However, the change of (v^) by the inclusion of the debris stars with f — 5 % 
amounts to only —19 (—11) km s^^ for q = 0.55 (0.7), which are still insufficient for explaining 
the reported (f^) = —35 ~ —65 km s~^. Also, if we extend the selection of the stars at higher 
z or instead consider model 1, the changes of (f,^) become smaller than the above mentioned 
values, because there are no debris stars in our current model. Thus, it is safe to conclude 
that the debris stars contribute only in part to a reported retrograde motion at the NGP. 

3.4. Local halo kinematics and microlensing towards LMC 

We select the nearby debris and halo stars at 6.5 < R < 9.5 kpc, z < A kpc, and D < 4 
kpc (as was drawn by CB), convolve the velocities with a Gaussian error distribution of 
1 cr = 30 km s"^, and compare with the corresponding stars with [Fe/H]< —2 in BOO. It 
follows that the non-Gaussian feature in velocities is much weaker than that at the NGP: 
the change of (v^) for / = 5 % is only —9 (—10) km s~^ in model 1 (model 2). This is due 
to the smaller debris fraction near z — than at high \z\. 

The effects of the debris stars on the microlensing optical depth towards LMC, r, are 
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modest as well. Following the Gould (1999) prescription for r and investigating the debris 
within 10° x 10° centered at LMC, we arrive at r < 10~^/, thereby indicating that r provided 
by the debris stars is much smaller than the observed r of O(10~^). 

4. DISCUSSION 

We have demonstrated that our fiducial models of an orbiting dwarf galaxy that once 
contained cu Cen predict a sequence of tidal streams in retrograde rotation and their existence 
is imprinted in kinematics of nearby stars, especially in the direction against Galactic rotation 
(GWN) and at the NGP (K02), while local halo kinematics and microlensing towards LMC 
remain unchanged. The streams are mostly distributed inside the solar circle, as suggested 
from the current orbit of uj Cen (Dinescu 2002). In contrast to the Sgr dwarf galaxy having 
polar orbit, the orbit of u Gen's progenitor galaxy is largely affected by a non-spherical 
disk potential, where the orbital plane exhibits precession with respect to the Galactic Pole, 
causing self-crossing of tidal streams in the disk. The projection of the orbit perpendicular 
to the disk shows an 'X'-like feature, thereby leaving denser streams at high l^j than at low 
\z\ for a given radius. This explains the significance of the debris at the NGP compared to 
the solar neighborhood. 

Existing kinematic studies of Galactic stars to search for a signature of u; Gen's progen- 
itor galaxy are yet confined to nearby stars, where the significance of the debris streams is 
modest, as shown here. Searches of stars inside the solar circle are more encouraging (Fig. 
1), in particular in the directions of / ~ 320° and I ~ 50°, where we expect the presence of 
high- velocity streams at Vios = 200 ~ 300 km s~^ and —400 ~ —300 km s~^, respectively. 
Future radial velocity surveys of these fields including the sample of the Sloan Digital Sky 
Survey or planned RAdial Velocity Experiment are worth exploring in this context. Also, 
detailed abundance studies of candidate stream stars will be intriguing, because such stars 
may exhibit different abundance patterns from field halo stars, as found in dwarf galaxies 
(Shetrone, Gote, & Sargent 2001). 

A yet unsettled issue is the origin of a progenitor satellite orbiting inside the solar 
circle, because dynamical friction alone from the present-day smooth Galactic components 
is insufficient for shrinking the orbit if it was born at a large distance (say, ~ 50 kpc) from 
Galactic center (Zhao 2002). One of the possibilities to preclude it may be that the merging 
of a satellite occurred while the Milky Way was still in the process of halo formation via 
hierarchical merging of several subgalactic clumps; successive gravitational interaction among 
clumps may help reducing the orbital angular momentum of a progenitor efficiently. Also, a 
progenitor may have formed in the vicinity of the proto- Galaxy, where the environment of a 
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strong tidal field promotes the formation of a compact M32-like galaxy (Burkert 1994) and 
its high density affords the survival chances until the epoch of the Galactic disk formation. 
Then, if a progenitor retained gas, growing Galactic tides induce the infall of gas into the 
progenitor center and trigger the formation of a globular cluster there (Bekki & Chiba 2002) 
similar to oj Cen. The story is yet speculative but worth pursuing based on sophisticated 
numerical codes for Galaxy formation. 

We are grateful to the anonymous referee for helpful comments. M.C. thanks Kenji 
Bekki and Tim Beers for useful discussions. 
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Fig. 1. — Upper (middle) panel shows the spatial distribution of the tidally disrupted debris 
for model 1 (model 2) after the 1.37 (1.86) Gyr orbital excursion of u Gen's progenitor 
galaxy. Lower panel shows the orbit of the galaxy center for model 1 (dotted hne) and 2 
(solid line). The plots are projected onto three orthogonal planes, where the Sun is located 
at x = — 8 kpc and xy corresponds to the disk plane. A frame measures 15 kpc on a side of 
each panel. The current position of lo Gen is at (x, y, z) — (—4.8, —4.0, 1.4) kpc. 



Fig. 2. — Velocity distribution of the debris particles in cylindrical coordinates, (dotted), 
(solid), Vz (dashed), for model 1 (a) and 2 (b). 



Fig. 3. — (a) Distribution of the heliocentric radial velocities in the direction of the GWN 
survey, for the debris stars of model 2 and randomly generated halo stars with g = 0.7. We 
select the stars at 1 < D < 5 kpc and 260° < I < 280° in the fields of 30° < |6| < 50° (with 
f — 1 % and 3 % for dotted and sohd histograms, respectively) and 20° < |6| < 40° (with 
f — 5%: dashed histogram), (b) The same as (a) but incorporating the randomly generated 
disk stars as well for 30° < |6| < 50° with / = 5 % (solid histogram). Dotted histogram 
denotes the contribution from the metal-poor halo alone. 
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